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Abstract 

Background: Targeted therapy against the Epidermal Growth Factor Receptor (EGFR) is among the most promising 
molecular therapeutics for Head and Neck Squamous Cell Carcinoma (HNSCC). However, drug resistance limits the 
clinical efficacy of anti-EGFR monoclonal antibodies and no predictive biomarker has entered the clinic yet. 

Methods: A retrospective clinical study was performed utilizing pathological specimens from 52 newly diagnosed 
HNSCC patients. These patients were screened for mutations in EGFR and KRAS. Tyrosine kinase mutations in EGFR 
and KRAS mutations were evaluated by high resolution melting analysis (HRMA), whereas EGFRvlll was determined 
using one-step real-time PCR. Finally, patient samples were screened for HPV-DNA by GP5+/6+ PCR. Survival analysis 
was performed using Kaplan-Meier analysis and significance was calculated using log-rank statistic. 

Results: In our study population no EGFRvlll mutations were present. However, two silent mutations were found; 
T785T in exon 20 and R836R in exon 21 of the EGFR gene. Additionally, HRMA revealed an abnormal KRAS melting 
pattern in 7.0% of the samples. However, the KRAS StripAssay could confirm only one sample with a G12S mutation 
and none of these samples could be confirmed by direct sequencing. HPV DNA was present in 3/25 larynx and 9/27 
oropharynx tumors. 

Conclusion: The low rate of EGFR and KRAS mutations in this Belgian HNSCC population suggests that these genes will 
probably not play a major role in predicting response to anti-EGFR therapy in HNSCC. Hence, other predictive markers 
need to be discovered in order to optimize EGFR targeting therapy. 

Keywords: Head neck squamous cell carcinoma, EGFR inhibitors, Cetuximab, Resistance, KRAS mutations, EGFRvlll 
mutations, EGFR tyrosine kinase mutations, HPV 



Background 

Head and neck squamous cell carcinoma (HNSCC) is a 
heterogeneous and complex disease, having a severe im- 
pact on the quality of life of patients and survivors [1]. 
At time of diagnosis, 60 to 70% of patients present with 
advanced disease, affecting the survival of these patients 
negatively [2]. In up to 90% of HNSCC tumors, in- 
creased expression of the epidermal growth factor recep- 
tor (EGFR) is observed [3-5]. Increased EGFR expression 
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is associated with advanced stage of the disease, poor sur- 
vival and resistance to chemo- and radiotherapy [3,6]. In 
only 50% of HNSCC patients, the current conventional 
treatment strategies, including surgery, chemotherapy and 
radiation, are effective, underscoring the need for new ap- 
proaches to treat this malignancy [7,8]. Moreover, the 
existing cytotoxic therapies are non-selective and associ- 
ated with considerable toxicity. Therefore, personalized 
medicine using targeted therapies will hopefully achieve 
the much- needed progress in HNSCC treatment. 

As EGFR signaling plays an important role in the de- 
velopment and progression of HNSCC, it is one of the 
most promising targets for molecular targeted therapies 
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[3-5]. The addition of cetuximab to platinum-fluorouracil 
chemotherapy, for instance, improved overall survival 
when given as first-line treatment in patients with recur- 
rent or metastatic HNSCC [9]. However, resistance to 
EGFR inhibitors presents a therapeutic challenge, since in- 
trinsic as well as acquired resistance after initial response 
has been described frequently [10]. Therefore, it is import- 
ant to understand the mechanism of intrinsic resistance, 
as a further selection of patients who might benefit from 
the addition of cetuximab will increase its beneficial effect 
and most likely its cost-effectiveness. 

In other forms of cancer where anti-EGFR therapy is 
used, molecular profiling of the tumor has become es- 
sential for customized medical treatment decision. EGFR 
and KRAS mutations strongly influence the benefit of 
treatment with anti-EGFR therapies, such as panitumu- 
mab and cetuximab used for treating metastatic colorec- 
tal cancer (mCRC) and tyrosine kinase inhibitors (TKIs) 
used for advanced non-small cell lung cancer (NSCLC) 
[11-14]. Although several reports have described EGFR 
and KRAS mutations in HNSCC, the mutation frequency 
can differ among ethnic groups. In the present study we 
aimed to correct for the heterogeneous nature of this 
disease by including only HNSCC patients eligible for 
cetuximab therapy and suffering from two common 
subtypes of HNSCC (oropharynx and larynx). More- 
over, since infection with human papillomavirus (HPV) 
is a driver of HNSCC carcinogenesis, and HPV-positive 
cancers form a distinct group with better clinical out- 
come within HNSCC [15,16], we determined the pres- 
ence of HPV infection as well. 

Methods 

Patient samples 

Formalin fixed, paraffin-embedded (FFPE) tissues from 
52 newly diagnosed HNSCC patient samples (27 oropha- 
ryngeal and 25 laryngeal SCC) were retrieved from the 
department of Pathology, Antwerp University Hospital 
(Edegem, Belgium). The tissues were obtained from bi- 
opsies or resections performed between 2002 and 2010. 
At that time, all patients agreed that residual material 
from resections was used for studies (consented at the 
time of hospitalization); no recurrent approval from the 
local ethic committee was required. As EGFR targeted 
therapies are approved in patients with locally advanced 
or recurrent/metastatic disease, our patients were se- 
lected for stage III or higher to most accurately resemble 
this population. 

Screening of EGFR tyrosine kinase and KRAS mutations 

Prior to DNA extraction, tumor tissue was enriched by 
manual macrodissection and isolated using the QIAamp® 
DNA FFPE Tissue isolation kit (Qiagen, Venlo, The 
Netherlands). Concentration and purity were defined 



using the Nanodrop ND-1000 spectrophotometer (Isogen, 
Sint-Pieters-Leeuw, Belgium). 

Genomic DNA from patients was used for mutation ana- 
lysis in the KRAS gene and EGFR tyrosine kinase domain 
by high resolution melting analysis (HRMA). Screening of 
mutations in exon 2 (codon 12 and 13) of the KRAS gene 
was performed as described previously [17]. Patient sam- 
ples with an abnormal melting pattern were subjected to 
KRAS StripAssay® (ViennaLab Diagnostics GmbH, Vienna, 
Austria), according to the manufacturer s recommendation. 
This assay screens for 10 KRAS mutations in codon 12 and 
13 of exon 2. 

Additional DNA was isolated from patients suspected of 
KRAS mutations, amplified by PCR and sequenced using 
Sanger sequencing. Primers for the 178 -bp amplicon of 
exon 2 were 5'-GTAAAACGACGGCCAGGTGTGACA 
TGTTCTAATATAG-3' (forward) and 5'-TTGGATCA 
TATTCGTCCACAA-3' (reverse). PCR reaction was per- 
formed in 25 (il reaction containing lx buffer, 1 mM 
MgCl 2 , 0.5 (iM of each primer, 0.25 mM dNTPs, 1 \A Taq 
polymerase and 5 [A genomic DNA (10 ng/|il). PCR cyc- 
ling was performed on a Px2 Thermal Cycler (Thermo 
Electron Corporation) and was run as follows: 15 minutes 
at 95°C and 35 cycles of 60 seconds at 94°C, 60 seconds at 
55°C and 60 seconds at 72°C, followed by 10 minutes at 
72°C. This PCR product was purified using ExoSAP-IT 
(Affymetrix, Cleveland, OH, USA) and for validation, this 
product was used as template for direct sequencing 
with the Big Dye Terminator vl.l kit (Applied Biosys- 
tems, Foster City, CA, USA) using M13tag primers 
(Eurogentec, Seraing, Belgium). The reaction mixture 
consisted of l.lx sequencing buffer, 0.2 \A Big Dye mix, 
625 nM M13tag primer 5'-GTAAAACGACGGCCAG- 
3' and 1 \A of purified template in a total volume of 
4 |lx1. The reaction was run on a Rapid Cycler Instrument 
2 (Idoha Technology Inc., Salt Lake City, UT, USA) ac- 
cording to the following protocol: initial denaturation 
at 95°C for 30 sec and 25 cycles in the following se- 
quence: 96°C for 10 sec, 50°C for 5 sec and 60°C for 
2 min. The sequencing reactions were run on a 3130 
XL Genetic Analyzer (Applied Biosystems). Sequencing 
data were analyzed using SeqScanner software vl.O 
(Applied Biosystems). 

Primers and reaction conditions for detection of EGFR 
tyrosine kinase mutations in exon 19, 20 and 21 were 
adapted from Heideman et al. [18]. Amplification and gen- 
eration of melting curves were performed on a Lightcycler® 
480 (Roche Diagnostics Gmbh, Mannheim, Germany) and 
analyzed using GeneScanning software (Roche). HRMA 
products with a deviating EGFR melting pattern were 
directly purified using ExoSAP-IT (Affymetrix) and the 
purified PCR product was used as template for direct 
Sanger sequencing (Big Dye Terminator vl.l kit (Applied 
Biosystems,)) using M13tag primers, as described above. 
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EGFRvlll mutation detection using real-time PCR 

The determination of EGFRvlll mRNA expression was 
performed by semi-quantitative real-time PCR. Total 
RNA was obtained from FFPE tumor tissue with the 
High Pure RNA Paraffin Isolation kit (Roche) according 
to the manufacturer s recommendation. Prior to extrac- 
tion, manual macrodissection was performed to enrich 
the sample for tumor cells. 

After extraction, one-step real-time PCR was performed 
on the LightCycler 480 I instrument (Roche) using Taq- 
Man RNA amplification kit (Roche). Thereby, 2 PCR reac- 
tions were performed, one for the specific amplification of 
EGFRvlll mRNA using a primer set targeting EGFR exon 
1 and exon 8 (EGFRel/8) and the other one for the ampli- 
fication of total EGFR mRNA using a primer set targeting 
an unaffected region of the EGFR gene (exon 9 and 10, 
EGFRe9/10). The primer pairs are given in Table 1. For the 
detection of the EGFRvlll product, a fluorescently labeled 
oligonucleotide TaqMan probe (Eurogentec) was used. For 
the detection of total EGFR PCR product, YOPRO inter- 
calating dye (Invitrogen, Merelbeke, Belgium) was used. 
Concentrated samples were diluted to a maximum of 
20 ng/ul. The input varied between 18.5 ng and 100 ng. 
RNA template was added in a reaction mix consisting of 
300 nM forward and reverse primer each, 3 mM Mn 2+ , 
0.1 mM YoPro (EGFRe9/10 reaction) or 200 mM TaqMan 
probe (EGFRel/8 reaction), lx RNA mix (TaqMan RNA 
amplification kit, Roche) in a total volume of 25 ul. Sam- 
ples were run in duplicate. 

As positive control, mRNA derived from an EGFRvlll 
mutant glioblastoma cell line U87del (kind gift from dr. 
Furnari from the Ludwig Institute for Cancer Research, 
La Jolla, CA, USA) was used. The cycling conditions are 
shown in Table 2. Data were analyzed using the Lightcycler 
480 software release 1.5.0 software (Roche). 

HPV detection 

HPV detection by GP5+/6+ PCR was carried out on all 52 
HNSCC samples, as described previously [19]. Detection 
of PCR products was performed in an enzyme immuno- 
assay format as described by Jacobs et al. [20]. 



of HPV infection was assessed by survival analysis using 
the Kaplan-Meier method and survival curves were ana- 
lyzed using the log-rank test. The index date for survival 
time calculation was defined as the date of diagnostic 
confirmation for HNSCC. The days of observation (over- 
all survival time) were calculated from the index date to 
the date of last information/ death. For progression free 
survival time, the days of observation were calculated 
from the index date to the first date of progression, 
death or the date of last information. All analyses were 
conducted using SPSS (version 20, SPSS Inc., Brussels, 
Belgium). Significance for all statistics was recorded if 
p < 0.05. 

Results 

Patient characteristics 

Of the 52 HNSCC patients from whom tumor tissue could 
be obtained, most clinical data could be retrieved. Further 
details on these patients are summarized in Table 3. 

EGFR mutation analysis: EGFR tyrosine kinase and 
EGFRvlll mutation 

HRMA was used as a first-line screening method for EGFR 
tyrosine kinase mutations. DNA samples yielding a suffi- 
ciently high DNA concentration and good purity were 
screened for mutations in exon 19 to 21 of the EGFR gene. 
In total, 46 HNSCC patient samples were available for EGFR 
mutation screening. Failure to generate a melting curve was 
present in 7, 7, 5 and 10 samples for exon 19, 20p, 20d and 
21, respectively. HRMA could not detect any aberrant melt- 
ing curves in exon 19, whereas in exons 20 and 21 abnormal 
melting curves were detected. Therefore, these samples were 
sequenced in order to confirm and determine the mutation. 
Although we did not detect any missense mutations, two 
patients were found with a silent mutation; T785T (exon 
20p) and R836R (exon 21) (Figure 1A and IB). 

Furthermore, screening for EGFRvlll mutations was 
performed by one-step real-time reverse transcriptase 
PCR. Of the 52 HNSCC samples, 8 failed to amplify. In 
the remaining 44 HNSCC samples, no EGFRvlll muta- 
tions could be detected. 



Statistical analysis 

All experimental data were the combined result of at 
least two independent experiments. Prognostic relevance 



KRAS mutation 

Similar to the screening method of EGFR tyrosine kinase 
mutations, we screened for mutations in codon 12 and 



Table 1 Primers (Eurogentec) used for detection of EGFRvlll by real-time PCR 


Oligo name 


Oligo type 


MocB 7 


Sequence 


Mod3' 


EGFRel-F 


forward primer 




GAGTCGGGCTCTGGAGGAA 




EGFRe8-R 


reverse primer 




GGCCOTCGCAOTCTOCA 




EGFRel/8-TM 


TaqMan probe 


6-FAM 


AAAGGTAATOTGTGGTGACAGATCACGGCTC 


BHQ-1 


EGFRe9/10-F 


forward primer 




TCCTGCCGGTGGCAm 




EGFRe9-R 


reverse primer 




TGTGGATCCAGAGGAGGAGTATG 
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Table 2 Cycling condition for one-step RT-PCR 



Table 3 Clinical data from 52 HNSCC patients 



Stage 



Temperature Time 



# cycles 



Oropharynx Larynx Overall population 



Stage I preheating 

Stage II: reverse 
transcription 



95°C 
60°C 



Stage III: activation 
Stage IV: cycling 



95°C 
95°C 
60°C 

Stage V: denaturation 95°C 
Stage VI: melting 50°C-85°C 
Stage VII: cool 40°C 



1 min 1 
30 min 1 

10 min 1 
15 sec 45 
1 min 
1 min 

Continues measuring 
30 sec 



13 of the KRAS gene using HRMA. Of the 52 HNSCC 
samples, five could not be used for screening due to re- 
duced DNA concentration or insufficient purity of the 
DNA samples. Another 4 samples failed to generate a 
melting curve; hence no mutation status could be deter- 
mined. Of the remaining 43 HNSCC patients, three 
samples presented with an aberrant melting profile, con- 
sistent with the positive control (KRAS mutated cell 
line), see Figure 2. The other 37 samples were wild type 
for codon 12 and 13 of the KRAS gene, as determined 
by HRMA. 

For validation, these three samples were subjected to 
further analysis using the KRAS StripAssay®. These re- 
sults showed a very weak band for the G12S mutations 
(green curve in Figure 2) for one sample, one other sam- 
ple was wild type for 10 KRAS mutations (yellow curve 
in Figure 2) and one sample could not be evaluated due 
to failure of the test (blue curve in Figure 2). To further 
confirm these results, additional DNA was isolated for 
KRAS mutation screening by dideoxy sequencing for 
these three samples. Surprisingly, all samples showed the 
wild type GGTGGC sequence. 

HPV infection in oropharygeal and laryngeal cancer 

Screening of our HNSCC patient population identified 12 
out of 52 HPV positive HNSCC tumors (23.1%), nine tu- 
mors were located in the oropharynx (33.3%) and three in 
the larynx (12.0%). 

Statistical analysis did not show any significant survival 
benefit (overall survival and progression free survival) for 
patients with HPV-related HNSCC tumors, p = 0.125 and 
p = 0.374 respectively (Figure 3 A and B). 

Discussion 

As personalized cancer medicine - in particular drugs 
based on genetics - will become more and more import- 
ant in the future, molecular profiling of the tumor prior to 
treatment becomes essential. Unfortunately, no predictive 
biomarker for the promising anti-EGFR therapeutics in 
the treatment of HNSCC has entered the clinic yet. 



Patient characteristics 








Total no. of patients 


27 


25 


52 


Median age (years) 


62.75 


62.17 


62.58 


Sex 








Male 


21 


21 


42 


Female 


6 


4 


10 


Grade of differentiation 








Good 


4 


7 


11 


Moderate 


17 


10 


27 


Poor 


0 


1 


1 


Unknown 


6 


7 


13 


TNM 








T1 


2 


0 


2 


T2 


6 


0 


6 


T3 


5 


10 


15 


T4 


14 


15 


29 


NO 


5 


12 


14 


N1 


1 


5 


6 


N2 


17 


6 


23 


N3 


3 


2 


5 


Nx 


1 


0 


1 


MO 


21 


22 


43 


M1 


3 


2 


5 


Mx 


3 


1 


4 


Stage 

III 


3 


8 


11 


Iva 


19 


13 


32 


Ivb 


4 


2 


6 


Ivc 


1 


2 


3 


Smoking 








Yes 


9 


6 


15 


Ex 


4 


6 


10 


No 


1 


0 


1 


Unknown 


13 


13 


26 


Alcohol abuse 








Yes 


7 


7 


14 


Ex 


1 


1 


2 


No 


3 


2 


5 


Unknown 


16 


15 


31 



Pinpointing mechanisms of resistance to EGFR target- 
ing agents in HNSCC has proven difficult due to several 
factors. First, HNSCC is a heterogeneous disease, affect- 
ing several anatomical structures located within the head 
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Figure 1 Sequencing result of two HNSCC patients. A: EGFR exon 20p results in a T> C mutation at nucleotide 2354 and B: EGFR exon 21 
displays at nucleotide 2508 C > T. Both mutations result in a silent nucleotide polymorphism, T785T and R836R for exon 20 and exon 21 respectively. 



and neck region. Second, mutation frequencies might 
vary between different ethnic groups. And finally, stages 
of tumors might be characterized by certain mutations. 
Therefore, we tried to correct for these factors by in- 
cluding only Caucasian laryngeal and oropharyngeal 



tumors eligible for cetuximab therapy and evaluated 
whether the mutation frequencies are in line with the 
frequencies reported in the literature. Furthermore, we 
determined the incidence of HPV in these 52 HNSCC 
patients as well. 



B 



65 
60 
55 
50 
45 
~U 40 

30 
25 
20- 
15 
10 
5- 
0- 
-5-, 




79 80 81 

Temperature 



<v 
u 
c 
<u 
u 
(/> 

<V 

o 



30 



25 



-10-, 




79 80 

Temperature 



Figure 2 Derivative plot (A) and difference plot (B) of the normalized high resolution melting curves of the three HNSCC samples 
suspected of KRAS mutation, together with two cell lines; A549 (homozygote G12S, red) and SC263 (homozygote wild type, black). All 

samples show an abnormal derivative plot compared to the wild type SC263 plot. 



Boeckx et al. BMC Research Notes 2014, 7:337 
http://www.biomedcentral.eom/1756-0500/7/337 



Page 6 of 8 



B 



Survival Functions 



Survival Functions 




■f- 0-censored 
I 1 -censored 



overall_survival 



S00 1000 1S00 2000 

progression_free_survival 



Figure 3 Kaplan-Meier survival analysis for HPV infection in HNSCC for overall survival (A) and progression free survival (B). Significance 
was calculated using log-rank statistic (p = 0.1 25 and p = 0.374, respectively). 



Sustained EGFR signaling can be elicited at different 
levels in the signaling cascade, including at the level of 
the target itself. Compensatory mutations in the tyrosine 
kinase domain of EGFR or the constitutively active and 
truncated form of EGFR (i.e. EGFRvIII) are shown to 
predict response to anti-EGFR therapy in other tumor 
types [13,14,21,22]. We observed no EGFR tyrosine kin- 
ase missense mutations, although two silent mutations 
(T785T and R836R) were detected. These results indi- 
cate that EGFR tyrosine kinase mutations are uncom- 
mon in this Belgian HNSCC patient population. This is 
in concordance with literature where mutation frequen- 
cies vary between 0 and 8% [23-25]. On the other hand, 
Murray et al. detected an EGFR tyrosine kinase mutation 
frequency of 15.8%, albeit in a small sample population 
(n = 19) [5]. The EGFRvIII mutation has been reported 
previously in HNSCC by real-time PCR and immunohis- 
tochemistry. However, subsequent studies have provided 
conflicting evidence [11,16,26,27]. In these studies, the 
EGFRvIII mutation frequency in HNSCC varies between 
0% and 42%. In our study population, no EGFRvIII mu- 
tations could be detected using one-step real-time PCR 
on FFPE material. 

KRAS serves as a mediator between extracellular lig- 
and binding and intracellular transduction of signals 
from EGFR to the nucleus. Therefore, aberrant expres- 
sion of KRAS might result in sustained EGFR signaling 
despite inhibition. The KRAS mutation frequency in 
HNSCC observed in literature varies between 0 and 6% 
[23,26,28]. This study identified three out of 43 samples 
with an aberrant HRMA melting profile. Although these 
samples showed aberrant melting profiles, indicative of 



the presence of low level sequence variants, no muta- 
tions could be identified with Sanger sequencing and 
only one (G12S mutation) could be confirmed with the 
KRAS AssayStrip®. These conflicting results can be ex- 
plained by two main factors, which influence test results. 
First, genetic heterozygosity within the tumor might re- 
sult in mutation levels below the analytical sensitivity of 
dideoxy sequencing (10 - 30%) [17,29], whereas the de- 
tection limit of HRMA was previously reported at 3% 
[17] and the KRAS StripAssay® detect at least 1% mu- 
tated alleles. Second, the DNA used in this study was ex- 
tracted from FFPE sections and formalin fixation has 
adverse effects on DNA. Consequently, Taq polymerase 
errors can cause PCR artifacts during amplification 
[29,30]. Therefore, positive samples were repeated by an 
independent round of DNA isolation and amplification, to 
avoid false positive results. However, it was not possible to 
make a distinction between low-level mutations and se- 
quence artifacts by formalin fixation, as fresh frozen tissue 
from these patients was unfortunately not available. 

As HPV positive HNSCC tumors tend to have a more 
favorable outcome, we screened our 52 patients for HPV 
DNA. In this Belgian study population, 33.3% of the oro- 
pharyngeal cancer and 12.0% of the laryngeal tumors 
contained HPV DNA. Frequencies reported in literature 
can range from 0% to 100% [31-33], depending on sev- 
eral factors, including differences in ethno-geographic 
regions, combination of different subtypes of HNSCC 
and differences in the applied HPV detection methods 
[34-36]. Our study was not designed to confirm the as- 
sociation between HPV positive HNSCC tumors and im- 
proved survival, and consequently, did not take into 
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account other factors influencing prognosis. Therefore, 
no conclusions regarding HPV positive HNSCC tumors 
and overall survival can be made. Furthermore, the value 
of cetuximab in HPV-related HNSCC is unclear. Recent 
data challenge the role of cetuximab therapy in HPV posi- 
tive oral squamous cell carcinoma [37,38]. Therefore, the 
role of cetuximab in treatment of HPV-associated HNSCC 
needs to be explored in prospective clinical trials. 

Conclusion 

To conclude, this retrospective study detected two silent 
EGFR tyrosine kinase mutations and no EGFRvIII muta- 
tions. Due to discrepancy between HRMA and sequencing 
results, KRAS mutation frequency varies from 0.0% to 
7.0%. Although the mutation status of EGFR and KRAS are 
used as predicting biomarkers for EGFR targeted therapy 
in other tumor types, their (extremely) low prevalence in 
HNSCC will likely preclude a major role in helping to de- 
fine treatment options in these patients. Therefore, other 
genes are likely involved in resistance to anti-EGFR therapy 
in HNSCC and more in-depth-studies are needed to un- 
ravel these mechanisms. 
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